function F = solve_f_mu(mu_f, sigma_f, cost_entry_ii, n_f_ii, N_f)

F = logncdf(cost_entry_ii, mu_f, sigma_f) - n_f_ii/N_f;

